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ABSTRACT 

We propose a new method to calculate Faraday rotation measure maps from multi-frequency 
polarisation angle data. In order to solve the so called n7t-ambiguity problem which arises 
from the observationally ambiguity of the polarisation angle which is only determined up 
to additions of ±n7T, where n is an integer, we suggest using a global scheme. Instead of 
solving the nTT-ambiguity for each data point independently, our algorithm, which we chose 
to call Pacman (Polarisation Angle Correcting rotation Measure ANalysis), solves the nn- 
ambiguity for a high signal-to-noise region "democratically" and uses this information to 
assist computations in adjacent low signal-to-noise areas. 

Key words: Intergalactic medium - Galaxies: cluster: general 



1 INTRODUCTION 



Magnetic fields are a common phenomenon on different astrophys- 
. ical scales. One possibility to observe astrophysical magnetic fields 
and thus, gain insight into their structure, ordering scale and field 
strength is given by the Faraday rotation effect. It arises when- 
• ever linearly polarised radio emission passes through a magnetised 
' medium. The polarisation angle ip of the radiation will be rotated 
, due to the interaction between the magnetic field component par- 
allel to the propagation direction of radiation and the radio emis- 
sion itself. The resulting change in polarisation angle is propor- 
tional to the wavelength squared which can be expressed by 
ip = RM A'^ -I- ip^ , where the proportionality constant RM in this 
relation is called Faraday rotation measure and the intrinsic 
polarisation angle of the emission at the radio source. 

The latter two quantities, RM and <^", can be measured 
through multi-frequency polarisation observations. The RM can be 
described as the line of sight integral over the electron density and 
the magnetic field component parallel to the line of sight. There- 
fore, the RM contains valuable information on the magnetic field 
in the foreground of the polarised radio source. The intrinsic polar- 
isation angle i^" and thus, the intrinsic magnetic field direction at 
the radio source gives insight into the magnetic field configuration 
at the source. 

The statistical analysis of RM measurements in terms of cor- 
relation fu nctions and equivalently power spectr a as developed by 
IVoet & En Blin (2003) and EnBlin & Voat; <2003h requires that the 
RMs are unambiguously determined. Thus, any ambiguous RM can 
lead to misinterpretation of the data investigated. For the calcula- 
tion of RM and ip'^ using the relationship p> = RM }? + pP men- 
tioned above, a least squares fit is normally applied to the polari- 
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sation angle data. Since the measured polarisation angle p> is con- 
strained only to values between and n leaving the freedom of ad- 
ditions of ± nn, where n is an integer, the determination of RM and 
^"is ambiguous, causing the so called riTt-ambiguity. Therefore, a 
least squares fit has to be applied to all possible nn-combinations 
of the polarisation angle data at each data point of the polarised ra- 
dio source while searching for the nyt-combination for which is 
minimal. 

In principle, can be decr eased to infin itely small num- 
bersby increasing RM substantiallv. lVallee & Kronbe rg ( 1975) and 
iHave 3 <1975h suggested to avoid this problem by introducing an 
artificial upper limit for |RM| < RMmax- Since this is a biased 
approach, Ruzmaikin & Sokolo proposed to assume that 

no riTt-ambiguity exists between the measurements of two closely 
spaced wavelengths taken from a whole wavelength data set. The 
standard error of the polarisation measurements is then used to con- 
strain the possible riTt-combinations for the least squares fit for the 
subsequent observed frequencies. We refer to these methods as the 
"standard fit" algorithms, as they are currently the most widely used 
methods. However, these methods might still give multiple accept- 
able solutions for data with low signal-to-noise, requiring the solu- 
tion to be flagged and all the information carried by these data is 
lost. Furthermore, it still can happen that the algorithm chooses a 
wrong RM and imprints spurious artefacts on the RM and the pp 
maps. 

Recently , a com pletely different approach was proposed by 
ISarala & laini i200lh which takes the circular nature of the polar- 
isation angle into account. The authors apply a maximum likeli- 
hood method to spectral polarisation data. Although this approach 
is not biased towards any RM value, it is rath er designed for a large 
number of observed wavelengths. Similarlv. lde Bruvi] (1996) and 
Brentjens & de Bruvn (2004) propose a RM-synthesis via wide- 
band low-frequency polarimetry. However, typically the observa- 
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tions are only performed at three or four wavelengths especially for 
extended (diffuse) radio sources. 

Here, we propose a new approach for the unambiguous deter- 
mination of RM and <^". We assume that if regions exhibit small 
polarisation angle gradients between neighbouring pixels in all ob- 
served frequencies simultaneously, then these pixels can be consid- 
ered as connected. Note, that we calculate the gradient modulo 7t, 
which implies that polarisation angles of and n are regarded as 
having the same orientation and thus, the cyclic nature of polarisa- 
tion angles is reflected. Information about one pixel can be used for 
neighbouring ones and especially the solution to the n7t-ambiguity 
should be the same. 

In cases of small gradients, assuming continuity in polarisa- 
tion angles allows us to assign an absolute polarisation angle for 
each pixel with respect to a reference pixel within each observed 
frequency. This assignment process has to be done for each spa- 
tially independent patch of polarisation data separately, such as 
each side of a double-lobed radio source. The reference pixel is 
defined to have a unique absolute polarisation angle and the algo- 
rithm will start from this pixel to assign absolute polarisation angles 
with respect to the reference pixel while going from one pixel to its 
neighbours. Figuratively, the algorithm eats its way through the set 
of available data pixel. It might become clear now why we have 
chosen to call the algorithm Pacman^ (Polarisation Angle Correct- 
ing rotation Measure ANalysis). 

Pacman reduces the number of least squares fits in order to 
solve for the riTt-ambiguity. Preferably, the reference point is cho- 
sen to have a high signal-to-noise ratio so that in many areas, it 
is sufficient to solve the nyt-ambiguity only for a small number 
of neighbouring pixels simultaneously and to use this solution for 
all spatially connected pixels. Pixels with low signal-to-noise will 
profit from their neighbouring pixels allowing a reliable determina- 
tion of the RM and ip° . 

In Sect. |2| we describe in detail the idea and the implemen- 
tation of the Pacman algorithm. In Sect.|3| we test Pacman on ar- 
tificially generated RM maps and demonstrated its ability to solve 
the nTT-ambiguity properly. However, the application to observa- 
tional polarisation data and statistical characterisation of the result- 
ing maps will be presented in the second paper (Vogt et al.), to 
which we will refer as Paper II. 

2 THE NEW PACMAN ALGORITHM 
2.1 The Idea 

As described in the introduction, the Faraday rotation measure 
RMij at each point with map pixel coordinate {ij) of the source, 
is usually calculated by applying a least squares fit to measured 
polarisation angles tpij(k) observed at frequency k G 1.../ such 
that 

VrAk)=RM,,\l+if%, (1) 

where iflj is the intrinsic polarisation angle at the polarised source. 

Since every measured polarisation angle is observationally 
constrained only to a value between and n, one has to replace 
ifiijik) in the equation above by f>ij{k) — fij{k) ± nij{k)n, 
where nij{k) is an integer, leading to the so called n7t-ambiguity. 
Taking this into account, a least squares fit to calculate RM^j and 

^ The computer code for Pacman is publicly available at 
fhttp : / / dipastro . pd .astro . it / ~cosmo/Pacman| 



ifiij at each pixel has to be applied by allowing all possible combi- 
nations of riij {k)n while determining the Uij (k) for which the Xij 
is minimal. The presence of observational noise might cause a stan- 
dard least squares fit, as suggested by Vallee & Kro nbera 
Haves ( 1975) or Ruzmaikin & S okolofj ] 1 97Sl) . to choose a spuri- 
ous RM value especially for areas of low signal-to-noise. 

The idea of Pacman is to reduce the number of pixels for 
which the n7t-ambiguity has to be individually solved. This is done 
by splitting the solution of the n7t-problem into two problems, a 
local and a global one, 

nij{k) = hij{k) + n{k), (2) 

where hij{k) is the local solution, linking polarisation angles of 
neighbouring pixels within a frequency map, and n{k) is the global 
solution to the problem, linking polarisation angles of the different 
frequencies. The local part hij (k) is determined by construction of 
absolute polarisation angle maps for each frequency with respect 
to a high signal-to-noise reference pixel being defined to possess 
a unique polarisation angle. The term absolute polarisation angle 
is to be understood as a value determined relative to the reference 
pixel by adding ztrni to the measured polarisation angles in order 
to remove jumps of the order of n in the measured polarisation 
angle map of each observed frequency. The global nn ambiguity 
is solved for a high signal-to-noise area surrounding the reference 
pixel resulting in n{k). This is then also the solution of the global 
riTt-problem for all spatially connected points which are assigned 
absolute polarisation angles with respect to this reference area. 

The splitting of the problem in a local and a global 
one is possible if the real polarisation angle ^Pij{k) 
(tpij{k) = KMij Xf. + (fi'^j) is a smooth quantity which does 
not change more than ±7t/2 between neighbouring pixels. In 
Sect. 12.21 it is described how Pacman ensures that only pixels 
fulfilling this condition are used. 

The source might consists of several spatially independent ar- 
eas of polarisation, to which we will refer as patches in the follow- 
ing. The nTt-ambiguity has to be solved for each of these patches 
separately, requiring separate reference pixels to be defined for the 
construction of absolute polarisation angles. 

The advantage of Pacman is that the global n7t-ambiguity is 
solved only for pixels having the highest signal-to-noise ratios. 
Therefore, noisier pixels which are situated at the margin of the 
source profit from an already defined global solution n{k) to the 
TiTt-ambiguity, making a reliable determination of RMij and the 
intrinsic polarisation angle (p°j for these pixels possible. 



2.2 The Basic Algorithm 

The algorithm starts with reading the maps of polarisation angles 
and of its errors, processes them, and when finished, it saves the 
different calculated maps. A flow chart of our Pacman algorithm 
is exhibited in Fig. Q which shows schematically the procedure 
the algorithm follows in order to determine the solution to the nn- 
ambiguity and to calculate the various maps. 

After loading of the various polarisation data for the different 
frequencies, the algorithm calculates a quality measure map. Dif- 
ferent actions are involved in the processing of the data. One of 
the processes is building a patch which involves the construction 
of absolute polarisation angle maps. First, the reference pixel for a 
patch of polarisation is determined from the quality measure map. 
This step is followed by solving the local rzTt-ambiguity in each 
frequency map for pixels being spatially connected to the defined 
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Figure 1. A flow chart representing the individual steps involved in the 
calculation of RM and ip" maps as performed by our algorithm Pacman. 



patch reference pixel and having small gradients in polarisation an- 
gle to the neighbouring pixels at all frequencies. This process is 
mainly governed and stopped by quality requirements. 

When the patch building process is stopped then the algorithm 
solves the global nyt-ambiguity across the frequency maps individ- 
ually for a selected set of superior quality pixels (voters) within the 
patch. The majority of votes is used as the global patch TiTT-solution 
and the resulting RM and tfP values are stored. If there are still un- 
processed high quality pixel left, the algorithm starts again to build 



a new patch. Otherwise Pacman finishes and saves the resulting 
complete maps. 

The values for RM and i/p" are calculated following a weighted 
least squares fit expressed by 



Sij \^ 



A2 



(3) 



(4) 



where S,, = Vc^i.,, Vrj = ELi f^ji^) /al^^, 

^^ELi^iKr AV7 = ELi^i<^>.«K, and 
— A|/cr|. . . In these relations, a^^j is the standard error 

of the polarisation angle fijik) at the pixel coordinate {ij) at the 
fcth wavelength Afc. 

Since the performance of the algorithm is mainly governed 
by quality requirements, we need some quality measure in order 
to rank the map pixels. One good candidate is the expected uncer- 
tainty aff^ of any RM value obtained in an error weighted least 
squares fit which is calculated by 



V S,, A* - A2' 

where the terms are defined as above. 

The uncertainty cr,^^ of the RM value is used to assign each 
pixel an initial quality qtj — since cr^^ accounts for the 

statistic used to determine the RM maps. Hence, small values of qij 
indicate high quality pixels. The quality is then used to determine 
the way Pacman goes through the data which will be preferably 
from high to low quality data pixels. 

For the construction of absolute polarisation angle maps for 
each frequency, Pacman starts at the best quality pixel having the 
smallest value qij, which is defined to be the reference pixel. For 
this first point, the measured polarisation angle ipij{k) is defined 
to possess a unique polarisation angle value (pijik) = 'Pij{k) for 
each observed frequency. It is important to note that the reference 
pixel is the same for all frequencies. 

Then, Pacman compiles what we call the border list (BL). It 
contains pixels being direct neighbours to points which have been 
already assigned an absolute polarisation angle. In the following, 
the set of adjacent pixels (i ± Ij ± 1) to the pixel (ij) will be 
referred to as direct neighbours (i'j'). Beginning with the direct 
neighbours of the reference pixels the border list is continuously 
updated during the progression of the algorithm. A pixel can be 
rejected from the border list if the standard error afc. of the polar- 
isation angle for any frequency at this pixel exceeds a certain limit 
(j™^" which is set at the beginning of the calculation. However, this 
requirement can be relaxed (see Sect. l2.7t . 

Having defined the reference pixel, Pacman assigns absolute 
polarisation angles to pixels within the border list always starting 
with the pixel having the best quality, i.e. the lowest qij value. For 
this pixel, the algorithm solves the local jiTT-ambiguity with respect 
to hij (k) £ Z hy minimising the expression 



^ [{ifiijik) ± h,j(k)n) - ipi,j'{k)f 



(6) 



{i'j'}SI3, 



where Dij is the set of all direct neighbours (i' j') to the pixel [ij) 
which have already been assigned an absolute polarisation angle. 
The resulting hij (k) value determines the absolute polarisation an- 
gle ^ij(k) — ipij{k) ± hij{k)n which has the smallest differ- 
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ence to the already defined absolute polarisation angles of adja- 
cent pixels ifiiiji (k). For each pixel, this is done at each frequency 
k € 1..../ independently but simultaneously. Thus, Pacinan goes 
the same way through the data in each frequency. Note that one can 
introduce a value a^^x which causes Pacman to reject the pixel if 
a exceeds (T^ax (see Sect. l2.6t . Additionally the border list is up- 
dated to include direct neighbours to the recently processed pixel 
{ij) which have not yet been assigned an absolute polarisation an- 
gle. 

Pacman repeats this process for the best remaining pixel in 
the border list which has the lowest qij and so on, until the whole 
patch consists of spatially connected pixels with assigned absolute 
polarisation angles and no acceptable neighbouring pixels remain 
in the border list. 

At this stage, Pacman solves the global nyt-ambiguity n{k) 
by applying a standard least squares fit to a set of best con- 
strained pixels within the patch (i.e. the surrounding area of 
the reference pixel). Pacman solves for each of these best con- 
strained pixels the global nTt-am biguity independentl y by either 
using the method oflHayesUl975l) :IVallee & KronberJ il975h or of 
iRuzmaikin & Sokolo fd <1979h minimising the for all possible 
combinations of nn. It searches then for the solution n(fc), which 
has been derived for the majority of pixels. This democratically 
determined solution n(fc) is used to perform a least squares fit of 
the other pixels defined within this patch. Note that in an extreme 
case, only the reference pixel could be used to solve the global nn- 
ambiguity. 

If there are still good quality pixels left in the map which were 
not connected to the first patch of good quality data, the algorithm 
begins a new patch of absolute polarisation angles, starting again 
with the best remaining pixel with smallest qij. A new patch of ab- 
solute polarisation angles is constructed. As a result, the source will 
be divided in a set of spatially independent polarisation patches, 
and the global n7t-ambiguity is solved once for each patch sepa- 
rately. A qfj^"^^^ can be introduced to prevent Pacman from starting 
new patches if the pixels remaining have quality values above that 
threshold value whose value is chosen in the beginning. 

2.3 Improving Quality 

The solution of the local n7t-ambiguity for pixel (ij) becomes more 
reliable as more direct neighbours of these pixels have their abso- 
lute polarisation angle ipiiji defined. This is taken into account by 
modifying the quality qij of the pixels added to the border list to 

where noij is the number of already defined neighbours (i'j'), and 
a and /3 are free parameters. This ensures that pixels having more 
already defined neighbours are considered before others although 
the original data point might have a smaller signal quality, i.e a 
higher qij. Values between and 1 for the free parameter P are 
reasonable and yield good results. We used a = (3 — 1 and ob- 
served that Pacman goes through the acceptable data points of the 
polarisation patches in a uniform manner. 

2.4 Restricting Gradients 

The algorithm might be faced with a situation where a high qual- 
ity pixel can be influenced by a very poor quality pixel, for ex- 



ample when the polarised radio source consists of two radio lobes 
each having a good signal-to-noise ratio which are connected by a 
bridge containing only low quality pixels. The Pacman algorithm 
would start by defining absolute polarisation angles from one of 
the two lobes, eventually reaching the bridge of low signal-to-noise 
and entering the second lobe from there. In such cases, it might 
happen that within the area of low quality data pixels a distinct 
determination of absolute polarisation angles by solving the local 
riTt-ambiguity using Eq. j5j is no longer possible. The entire second 
lobe would then suffer, and wrong solutions would be introduced. 

In order to avoid such situations, the algorithm is restricted 
to accept in the border list only neighbours which have a lower 
quality than the one under consideration. This forces the algorithm 
to go always from high to low quality pixels, leading to an artificial 
splitting of connected regions in the map into different patches. 

However, such a strict rule would lead to heavy fragmentation 
and is not preferable. Therefore, we introduced a parameter g to 
relax this strict limitation such that a new pixel is only accepted in 
the border list when the relation 

Qi'j' > 9* Qij (8) 

is fulfilled. This relation is always accounted for when adding a 
new pixel to the border list. We found that g between 1.1 and 1.5 is 
a good choice. 



2.5 Topological Defects 

The algorithm can also be faced with situations which we call topo- 
logical defects. These defects can be understood by supposing a 
ring-like polarised structure for which all polarisation angles point 
towards the centre of the ring. Starting at any point of the ring to 
define absolute polarisation angles and following the ring struc- 
ture, a jump will appear on the border of the first to the last defined 
absolute polarisation angles after having performed a full circle. 
IScheuer et aljil977l) describe this problem which they encountered 
by their analysis of polarisation data. 

In observational data, topological defects are often more com- 
plex structured. The artificial jump introduced will cause difficul- 
ties in the solution of the local 7i7t-ambiguity because this proce- 
dure relies on already defined absolute polarisation angles of neigh- 
bouring pixels. When confronted with this situation, our algorithm 
divides the list of direct neighbours into sublists possessing similar 
polarisation angles. The sublist containing the most pixels is then 
used to assign the absolute polarisation angle by solving the local 
n7t-ambiguity for the pixel under consideration. 

The locations of the polarisation angle steps of topological de- 
fects are somewhat artificial since they depend on the actual path 
of the algorithm through the data. Since the algorithm processes 
all frequencies simultaneously, the steps of topological defects are 
at the same positions within all frequency maps, and therefore do 
not cause any further problems. However, if Pacman encounters a 
jump in the polarisation angle to all possible neighbour sublists in 
any of the frequencies under consideration, Pacman will reject this 
pixel and this pixel will not be considered for this patch but queued 
back for consideration for the next patches to be constructed. In 
our experience, these topological defects are rare events, but it is 
necessary to take them into account for the solution of the local 
riTt-ambiguity. 
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2.6 Spurious Points 

The intrinsic polarisation angle distribution might show strong gra- 
dients extending over a few pixels. This could lead to situations, 
where our algorithm can not solve the local riTt-ambiguity at all 
frequencies simultaneously. Therefore, Pacman refuses to assign 
an absolute value (pij to pixels when Eq. (4) yields a value cr^ 
above a certain threshold cr^^^, which can be set at the beginning 
of the calculation. 

The pixels in the regions where this might occur most often 
have a low signal-to-noise ratio. In our experience, such situations 
always occur at strongly depolarised areas, leading to blanked re- 
gions in the RM distribution. 

2.7 Multi-Frequency Fits 

The aim of any RM derivation algorithm should be to calculate 
RMs for an area as large as possible using as much information 
as is available. On the other hand, for radio observations the total 
radio intensity decreases with increasing frequency. This can lead 
to the problem that the area of acceptable polarisation data at a 
high frequency is much smaller than at a lower frequency. This is 
especially true for (diffuse) extended radio sources. Furthermore, 
the limit of allowed standard errors cr^'^^ of the polarisation angle 
might be exceeded for only one frequency leaving the values for 
the other frequencies still in the acceptable range. 

Pacman can take this into account and performs the RM fit by 
omitting the polarisation angles at frequencies which do not meet 
the quality requirements. In order to do that, Pacman uses the stan- 
dard errors crfc . of polarisation angles to define independently for 
each frequency, areas in which the RM map is permitted to be pro- 
duced. An additional parameter fcmin = / is introduced which de- 
scribes the minimum number of frequencies allowing the freedom 
to use any combination of the minimum of frequencies. 

In C3.SCS of fcmin < /, Pacman will start to determine the solu- 
tion to global and local nyt-ambiguity only for the pixels fulfilling 
the quality requirements at all / frequencies. After finishing that, 
the algorithm proceeds to include pixels satisfying the quality crite- 
ria in less than / frequencies. For these pixels, the same patch build- 
ing procedure as described in Sect. 12.21 applies with some modifi- 
cation in order to prevent the final maps from heavy fragmentation. 
The best quality pixel among the remaining pixel is picked but be- 
fore starting a new patch, Pacman tests if the pixel under consider- 
ation adjoins a patch which has already been processed before. If 
the pixel adjoins such a patch, Pacman tries to solve the local nn- 
ambiguity following Eq.j6j and applies the patch solution of the 
global nTt-ambiguity to the pixel. If the pixel is neither adjoined 
to an already processed patch nor the local riTt-ambiguity solvable 
(i.e. cr^ exceeds cr,'^ax) then a new patch is initialised by this pixel. 

Another possibility is to force Pacman to use certain frequen- 
cies in all circumstances. Thus, if the quality requirements are not 
fulfilled for these particular frequencies, there will be no RM value 
determined for the pixel under consideration. This has the advan- 
tage that one can use relatively close frequencies as a basis and 
then include other frequencies at points when a reliable polarisa- 
tion signal is detected. The advantage of this technique is discussed 
in Paper II. 

2.8 Additional Information 

Apart from the resulting RM and (p^ maps, Pacman provides sets of 
additional information about the data in order to estimate the relia- 



bility of the results. A patch map which contains all patches used is 
one example. Such a map is very useful, especially if one requires 
a minimal number of pixels in a patch in order to accept any calcu- 
lated RM values from a particular patch. A map which includes re- 
jected and thus flagged pixels, can also be obtained. Probably more 
important are the final y^- and a^^- maps which are also provided 
by Pacman. This information allows one to understand the reliabil- 
ity of the RM maps obtained and can be used for further evaluation 
and analysis of the RM maps. 



3 TESTING THE PACMAN ALGORITHM 

In order to demonstrate the ability of Pacman to solve the titt- 
ambiguity, we use artificially generated maps. To g enerate them, 
we start from a polarisation data set from Abell 2255 il Govoni et alj 
2002), which was kindly provided by Federica Govoni. We calcu- 
lated the RM and (/jq maps from this data set and assume that these 
two maps consist of exact values. We generated then polarisation 
maps at four different frequencies which would exactly result from 
the RM and tpo maps. As error maps, we generated Gaussian devi- 
ates for each frequency which were multiplied by the original error 
maps provided by Federica Govoni. The so generated error maps 
were added to the generated frequency maps in order to provide 
realistic mock observations. 

Pacman and the standard fit algorithm were used to calculate 
the corresponding 'observed' RM and <^o maps which then were 
compared pixel by pixel to the initial exact maps. The result of this 
comparison is shown in Fig.|2| On the right panel, the comparison 
between the RMpacman values of the Pacman map and the val- 
ues RMsim of the initial map is shown. The scatter in the data is 
due to the noise which was added to the frequency maps. On the 
left panel of Fig.|2| the pixel by pixel comparison between the val- 
ues of the standard fit RMstdflt map and the initial RMsim map is 
shown. Again the scatter is due to the added noise. However, one 
can clearly see the points at ± 1000 rad m~^ which deviate from 
the initial data and are due to the wrongly solved nvr-ambiguity in 
the case of the standard fit. Thus, this test demonstrates that Pacman 
yields reliable results to the nyr-ambiguity. 



4 CONCLUSIONS 

We have presented a new algorithm for the calculation of Faraday 
rotation maps from multi-frequency polarisation data sets. Unlike 
other methods, our algorithm uses global information and connects 
information about individual neighbouring pixels with one another. 
It assumes that if regions exhibit small polarisation angle gradients 
between neighbouring pixels in all observed frequencies simulta- 
neously, then these pixels can be considered as connected, and in- 
formation about one pixel can be used for neighbouring ones, and 
especially, the solution to the n7t-ambiguity should be the same. 
We like to stress, that this is a very weak assumption, and - like 
all other criteria used within the Pacman algorithm - only depends 
on the polarisation data at hand and the signal-to-noise ratio of the 
observations. 

Our Pacman algorithm is especially useful for the calculation 
of RM and maps of extended radio sources. Global algorithms 
as implemented in Pacman are preferable for the calculation of RM 
maps and needed if reliable RM values are desired from low signal- 
to-noise regions of the source. 

Our Pacman algorithm reduces n7t-artefacts in noisy regions 
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Figure 2. Comparison between the simulated RMgim values for an artificially generated set of polarisation angle maps of four frequencies and values derived 
by Pacman RMpacman and the standard fit algorithm RMstdfit- The scatter is due to noise which was added to the polarisation angle maps. Note the wrong 
solutions of the n7r-ambiguity results which were calculated by the standard fit indicated by the points around ± 1000 rad m"^. 



and makes the unambiguous determination of RM and iy3° in these 
regions possible. Pacman allows to use all available information 
obtained by the observation. Any statistical analysis of the RM 
maps will profit from these improvements. 

A more detailed description of the application of Pacman to 
observations and various statistical tests performed on resulting 
RM maps in comparison to the standard fit algorithms are given 
in Paper II. 
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